arXiv:1205.4402v2 [astro-ph.HE] 2Jul2012 


Mon. Not. R. Astron. Soc. 000, 1-12 (2012) Printed 3 July 2012 (MN MgX style file v2.2) 


Energy-dependent orbital modulation of X-rays and constraints on 
emission of the jet in Cyg X-3 

Andrzej A. Zdziarski, 1 * Chandreyee Maitra , 2,3 Adam Frankowski 1,4 , 

Gerald K. Skinner 5,6 and Ranjeev Misra 7 

1 Centrum Astronomiczne im. M. Kopernika, Bartycka 18, PL-00-716 Warszawa, Poland 

-Raman Research Institute, Sadashivnagar, Bangalore 560080, India 

'Joint Astronomy Programme, Indian Institute of Science, Bangalore 560012, India 

4 IS DC Data Centre for Astrophysics, University of Geneva, ch. d’Ecogia 16, CH-1290 Versoix, Switzerland 

5 Astroparticle Physics Laboratory, Code 661, CRESST and NASA/Goddard Space Flight Center, Greenbelt, MD 20771, USA 

6 Department of Astronomy, University of Maryland, College Park, MD 20742, USA 

1 Inter University Centre for Astronomy and Astrophysics, Pune University Campus, Pune 411007, India 


Accepted 2012 June 29. Received 2012 June 28; in original form 2012 May 20 


ABSTRACT 

We study orbital modulation of X-rays from Cyg X-3, using data from Swift , INTEGRAL and 
RXTE. Using the wealth of the presently available data and an improved averaging method, 
we obtain energy-dependent folded and averaged light curves with unprecedented accuracy. 
We find that above ~5 keV, the modulation depth decreases with the increasing energy, which 
is consistent with the modulation being caused by both bound-free absorption and Compton 
scattering in the stellar wind of the donor, with minima corresponding to the highest optical 
depth, which occurs around the superior conjunction. We find a decrease of the depth below 
~3 ke V, which appears to be due to re-emission of the absorbed continuum by the wind in soft 
X-ray lines. Based on the shape of the folded light curves, any X-ray contribution from the jet 
in Cyg X-3, which emits y-rays detected at energies >0.1 GeV in soft spectral states, is found 
to be minor up to ~ 1 00 keV. This implies the presence of a rather sharp low-energy break in 
the jet MeV-range spectrum. We also calculate phase -resolved RXTE X-ray spectra, and show 
the difference between the spectra corresponding to phases around the superior and inferior 
conjunctions can indeed be accounted for by a combined effect of bound-free absorption in 
an ionized medium and Compton scattering. 

Key words: radiation mechanisms: non-thermal - stars: individual: Cyg X-3 - stars: winds, 
outflows - X-rays: binaries. 


1 INTRODUCTION 

Cyg X-3 is a high-mass X-ray binary with a Wolf-Rayet (WR) 
donor (van Kerkwijk et al. 1992, 1996; van Kerkwijk 1993) and a 
very short orbital period of P - 4.8 h. Its distance is = 7-9 kpc 
(Ling, Zhang & Tang 2009; Dickey 1983; Predehl et al. 2000). In 
spite of its discovery already in 1966 (Giacconi et al. 1967), Cyg X- 
3 remains poorly understood. In particular, the nature of its compact 
object remains uncertain, due to the lack of reliable determination 
of the mass functions and inclination (see, e.g., Vilhu et al. 2009 
for a discussion). However, the presence of a black hole is favoured 
by considering the X-ray and radio emission and the bolometric 
luminosity (Hjalmarsdotter et al. 2008, 2009; Szostek & Zdziarski 
2008, hereafter SZ08, Szostek, Zdziarski & McCollough 2008, 
hereafter SZM08). Also, Zdziarski, Misra & Gierlinski (2010) 
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have shown that the differences in the form of the X-ray spectra 
of Cyg X-3 from those of confirmed black-hole binaries can be ac- 
counted for by Compton scattering in a cloud formed by the stellar 
wind from the companion. 

Cyg X-3 is a persistent X-ray source with a typical X-ray lu- 
minosity of Lx ~ 10 38 erg s -1 . Its X-ray spectra have been classified 
into five states by SZM08, who have also quantified their correla- 
tions with the radio emission. Its high-energy y - ray emission has 
been discovered by the Fermi Large Area Telescope (LAT) and by 
AGILE in the soft spectral states (Fermi LAT Collaboration 2009, 
hereafter FLC09; Tavani et al. 2009). Later detections by the LAT 
and AGILE are presented by Williams et al. (2011), Corbel et al. 
(2012) and Bulgarelli et al. (2012), hereafter B12. 

Cyg X-3 shows pronounced flux modulation on the 
4.8 h period, discovered in X-rays (Parsignault et al. 1972; 
Sanford & Hawkins 1972; Canizares et al. 1973) and in in- 
frared (Becklin et al. 1973). Both periodicities are most likely 
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related to the strong stellar wind from the WR star (e.g., 
Willingale, King & Pounds 1985; van Kerkwijk 1993). In particu- 
lar, the X-ray flux minima correspond to the maximum absorption 
in the line of sight around the superior conjunction (the compact 
object behind the WR star). The X-ray period has been found to 
increase (Manzo, Molteni & Robba 1978; Lamb, Dower & Fickle 
1979; Eisner et al. 1980; van der Klis & Bonnet-Bidaud 1981, 
1989; Kitamoto et al. 1987, 1995; Singh et al. 2002). The period 
increase of P/P - 1 x 10' 6 yr' 1 is likely to be due to the loss of an- 
gular momentum through the wind mass loss (Davidsen & Ostriker 
1974). 

In this work, we present a detailed study of the dependence 
of the form of the X-ray orbital modulation on energy, using data 
from Swift, INTEGRAL and Rossi X-ray Timing Explorer ( RXTE ). 
We also study the dependence of the modulation on the X-ray spec- 
tral state, as well as obtaining X-ray modulation profiles for inter- 
vals simultaneous with the detections of high-energy y-rays. The 
y- ray emission was found to be also strongly modulated on the or- 
bital period, but with the maxima and minima approximately corre- 
sponding to the X-ray minima and maxima, respectively (FLC09). 
This phase reversal has been interpreted as due to anisotropy of 
Compton scattering of the stellar emission by relativistic electrons 
in the jet (Dubus, Cerutti & Henri 2010). The jet emission certainly 
extends to lower energies, and, using the model of Dubus et al. 
(2010), the hard X-rays from the jet are also predicted to have the 
maximum around the superior conjunction (Zdziarski et al. 2012, 
hereafter Z 12), i.e., shifted with respect to the observed X-ray mod- 
ulation by a half of the period. A measurement of the energy- 
dependent orbital modulation up to hard X-rays simultaneous with 
the y- ray emission can then constrain the X-ray contribution of the 
jet. 


2 FOLDED LIGHT CURVES FROM X-RAY 
MONITORING 

2.1 The X-ray monitoring data 

We use X-ray monitoring data from the Swift Burst Alert Tele- 
scope (BAT; Barthelmy et al. 2005; Markwardt et al. 2005) in 
the form of a 14-195 keV 8-channel light curve created for 
this study. The typical exposure for a single data point is ~ 
10 2 3 s. The data contain 23783 points within 1807 days span- 
ning MJD 53355-55469. The channels are between energies of 
14, 20, 24, 35, 50, 75, 100, 150 and 195 keV. These data 
have been obtained using an analysis different from that of 
http://swift.gsfc.nasa.gov/docs/swift/results/transients, which gives 
BAT count rates in the single 15-50 keV channel. Fig. 1 shows 
a comparison of the daily-averaged 15-50 keV count rates with 
the corresponding 14-50 keV energy fluxes from our analy- 
sis. The energy fluxes have been obtained by converting the 
count rates of our data using scaling to the Crab spectrum (see 
Zdziarski, Pooley & Skinner 201 1). 

We also use both dwell and daily-averaged data from 
the All-Sky Monitor (ASM; Bradt, Rothschild & Swank 1993; 
Levine et al. 1996) on board RXTE. The dwell data contain 97996 
measurements within 5267 days spanning MJD 50087-55915, and 
the exposure of a single observation is < 90 s. The ASM has three 
channels at energies of 1.5-3 keV, 3-5 keV and 5-12 keV, for which 
the corresponding energy fluxes are obtained by scaling to the Crab 
spectrum, as above. 

Fig. 2 shows the relationship between the energy fluxes in 



Figure 1. The correlation between the daily averages of the BAT flux, F(14- 
50 keV) from our analysis vs. those of the 15-50 keV rate in the public BAT 
data, see Section 2.1. Only points with statistical significance > 3<r (for each 
of the quantities) are shown. The black cross shows the flux/rate of F(14-50 
keV) = 1.2 keV cnT 2 s' 1 , corresponding to 2?( 1 5—50 keV) = 0.028 cm' 2 
s' 1 , below which we define our soft spectral state. 
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Figure 2. The correlation between the daily-averaged energy fluxes in the 
3-5 keV and 14-50 keV ranges. The hard and soft states are defined here by 
F( 3-5 keV) < 0.58 keV cm' 2 s' 1 and F(14-50 keV) < 1.2 keV cm' 2 s' 1 , 
respectively. The boundaries are marked by the dotted and dashed lines. An 
intermediate state corresponds to both the 3-5 keV and 14-50 keV fluxes 
above the respective boundaries. Only points with statistical significance 
> 3<r (for each of the fluxes) are shown. 
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the 3-5 keV and 14-50 keV photon energy ranges. We see that 
at low soft X-ray fluxes, about F( 3-5 keV) < 0.58 keV cirT 2 s -1 
(marked by the dotted line), corresponding to the ASM 3-5 keV 
count rate < 3 s _1 , there is a clear soft/hard X-ray flux anticorre- 
lation. Note that a positive correlation between the 3-5 keV count 
rate and radio flux is seen below the same rate (SZM08). At higher 
3-5 keV fluxes but at high 14-50 keV fluxes (F(14— 50 keV) >1.2 
keV cut 2 s , marked by the dashed line), the hard X-ray flux 
changes within a narrow range independent of the soft X-ray flux. 
At lower 14-50 keV fluxes (below the dashed line), there is an ap- 
parent anticorrelation down to F( 14-50 keV) ~ 0.5 keV cuT 2 s _1 . 
Below it, the soft X-ray flux changes within a narrow range approx- 
imately independent of the hard X-ray flux. Such anticorrelations, 
expressing a spectral pivoting around ~ 10-15 keV, can also be in- 
ferred from figs. 1-2 of SZM08. Fig. 7b of SZM08 and fig. 3 of 
Hjalmarsdotter et al. (2008) show similar correlations with the 20- 
100 keV flux measured (though with a lower sensitivity than that of 
the BAT) by the Burst and Transient Source Experiment (BATSE) 
on board Compton Gamma Ray Observatory ( CGRO ). 

We have also studied MAXI (Matsuoka et al. 2009) data. 
However, the exposure of a single data point for those data is 1.5 h, 
which is ~ P/3 and much longer than any orbital bin we use. Thus, 
we do not use those data for studying orbital modulation in Cyg X- 
3. We note those data also provide us with a hard-state condition, of 
the 2-4 keV photon flux <0.2 cuT 2 s -1 , approximately equivalent 
to the ASM 3-5 keV count rate < 3 s _1 . 

We perfonned analyses including all of the ASM data and sim- 
ilarly with all of the BAT data, as well as some subsets of each. 
First, we use a hard-state criterion of the daily-averaged 3-5 keV 
ASM count rate < 3.0 s -1 , corresponding to P(3-5 keV) < 0.58 
keV cm -2 s -1 , to the left of the dotted line in Fig. 2. As the soft- 
state criterion, we use the daily-averaged 14-50 keV BAT flux of 
<1.2 keV cuT 2 s -1 , i.e., below the dashed line in Fig. 2. Our inter- 
mediate state corresponds to both fluxes being above these bound- 
aries, i.e., in the upper right comer of Fig. 2. The soft-state cri- 
terion corresponds to the intervals of MJD 53746-53896, 5391 3— 
54138, 54204-54252, 54445-54470, 54548-54642, 54748-54818, 
54980-55044, 55080-55103, 55322-55349 and 55584-55645 (de- 
termined using a 5-d running average in order to reduce the effect 
of fluctuations). For epochs after the end day of our BAT data, we 
have applied the corresponding criterion for the public 15-50 keV 
BAT data. 

We also use data from intervals corresponding to the detected 
high-energy y- ray emission. Currently, they are MJD 54566-54647 
(B12), 54750-54820 (FLC09), 54821-54850 (B12), 54990-55045 
(FLC09), 55324-55326 (B12), 55343-55345 (Williams et al. 
2011; B12) and 55586-55610, 55642-55644 (Corbel et al. 2012). 
All those periods correspond to the soft state as defined above, i.e., 
below the dashed line in Fig. 2, see FLC09 and B12. Thus, the ap- 
pearance of y-ray emission corresponds to a marked change of the 
form of the soft/hard X-ray correlation. Combining the results from 
Fermi LAT and AGILE, it appears likely that all low hard-flux in- 
tervals are also associated with some high-energy y-ray emission. 

2.2 Treatment of the light curves 

Cyg X-3 has an orbital period that is increasing with time. Thus, 
we use a quadratic ephemeris, 

T m = T 0 + P 0 m + cm 2 , c = P 0 P/2, P = P 0 +2cm, (1) 

where T m is the time of a 777-th occurrence of a zero orbital phase 
(presumably related to the superior conjunction) since the reference 


time, To, Po is the period at To, P is the period derivative, and P is 
the period at T„,. The ephemeris (in UTC) of Singh et al. (2002) is 

T 0 = 40949.392 [MJD], P 0 = 0.19968443 d, c = 5.75 x 10~ u d.(2) 

An updated ephemeris taking into account Suzaku observations (S. 
Kitamoto, private communication) is 

T 0 = 40949.3913 [MJD], P 0 = 0.19968451 d, c = 5.62x10-“ d,(3) 

in which case P/P =* (1.03 ± 0.02) x 10~ 6 yr _1 . 

We note that these ephemerides of Cyg X-3 use the template 
of van der Klis & Bonnet-Bidaud (1989), which is defined numeri- 
cally by their table 2. That template has the minimum at the phase 
— 0.96-0.97 instead of 1.0 (or, equivalently, phase 0.0). Thus, min- 
ima at phases < 1 obtained by us do not indicate a discrepancy 
with respect to previous results. We also correct the light curves for 
barycentric delays. 

We need to properly average the folded light curves, [F, ■(/■), 
<Tj(l;)\, where F t is the average count rate of the observation with 
the mid-time of tj, and <x,- is its measurement error. One issue to 
consider is the large dynamic range of the flux variability of Cyg 
X-3, which is almost two orders of magnitude (see Fig. 2). Thus, 
linearly-averaged modulation profiles are strongly biased by those 
at highest flux states. Using flux logarithms can, in principle, some- 
what alleviate this problem, as it strongly reduces the dynamic 
range used in summation. However, this requires removal of all 
negative flux measurements. For a weak signal, this introduces a 
strong bias against the minima of the modulation. We found that 
this effect strongly suppresses the depth of the modulation at ener- 
gies > 30 keV. 

On the other hand, orbital modulation acts on a signal with 
a local flux level, the level being governed by the source aperiodic 
variability. The periodic and aperiodic variabilities are largely inde- 
pendent of each other. Thus, in order to effectively reduce the large 
amplitude of the underlying aperiodic flux variability, we calculate 
a running average, averaging all observations within some ±A t of a 
given observation time, 

Fi = jYi F L (4) 

i 

where J is the number of data points satisfying the above condi- 
tion on j. Then, we renormalize both F, and cr, by dividing each by 
Fj (so (T,/F, remain unchanged). In this way, effects of the aperi- 
odic long-term variability on time scales > 2A t are removed. In the 
case of measurements equally spaced at a time interval of <K P/K, 
where K is the number of bins per period at which the folded light 
curve is calculated, using A t = P/2 would optimally remove the ef- 
fect of long-term aperiodic variability. However, our measurements 
do not satisfy this criterion, and are occasionally sparsely and un- 
evenly spaced. Then, the shorter the At, the worse the statistics on 
which the running average is based. In the extreme case of no other 
measurements within ±A t from a given observation time, F; would 
be reset to unity, and such a point would falsely reduce the actual 
orbital modulation. Thus, in our method, we reject points for which 
the running average is based on J < 10 points. We have tested this 
method for our data, and have found that At = 1 d (- 5 P) is ap- 
proximately the shortest interval for which the running average is 
estimated with good statistics and only a small fraction of points 
have J < 10. Namely, ( J ) ± rms, and the fraction of the rejected 
points are (51 ± 23, 0.02), (34 ± 16, 0.04), for the ASM and BAT 
data, respectively. We thus use the above value of At hereafter. 

Another issue here is related to aperiodic short time-scale vari- 
ability on time scales shorter than the length of an orbital bin, which 
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effect we also would like to remove. Therefore, we pre-average the 
light curves on real-time intervals with the length equal to the bin 
size, P/K, where K is the number of phase bins per orbit. We obtain 
local averages of the renormalized (see above) flux and its error, 
which we denote as F kI , a kh respectively, where k is the number 
of the phase bin and / is the number of the time bin (counted from 
the start time of a given data set) contributing to the k-th phase bin. 
This method was applied in studying periodic variability of Cyg 
X-l by Ibragimov, Zdziarski & Poutanen (2007). It also partly re- 
moves a bias on folded averages due to a non-uniform coverage, 
e.g., due to a large number of points during one time interval and a 
low number of points during another, with both contributing to the 
same phase bin, k. Without pre-averaging, the former would have 
much higher weight than the latter, even if all measurement were 
relatively accurate. The pre-averaging does not have a major effect 
in the case of Cyg X-3 due to its very short orbital period (result- 
ing in rather few points per a P/K time interval), as well as for 
ASM and BAT data, which have relatively unifonn coverage. Still, 
this is the statistically correct procedure, resulting, in particular, in 
correctly estimated uncertainties of the folded and averaged light 
curves. 

In final averaging, we use weights given by inverse squares of 
the uncertainties, 


(F k ) 




(5) 


where N k is the number of time bins contributing to the £-th phase 
bin. The error of the weighted average is estimated from the vari- 
ance, ((F 2 k ) - (F k ) 2 )/(N k - 1), 


<«*> = 


N k - 1 


( x'Nk tt ’2 i—l 

Z/=i 

y Nt Her 1 

, Zj/ = i 


■ < F k y 


( 6 ) 


We fit the obtained folded/averaged light curves by a sim- 
ple model of absorption/scattering in a spherically symmetric stel- 
lar wind. A commonly-used wind velocity profile from a massive 
star (e.g., Lamers, Cerruti-Sola & Perinotto 1987) and the resulting 
electron density are, respectively, 


/ R, / -M 

v(r)^v„ 1 , n e (r )=- 7 2 " 7 v ( 7 > 

where r is the distance from the centre of the donor of the radius R„ 
v m is the terminal velocity, /? parametrizes the wind acceleration, M 
is the total mass-loss rate, m p is the proton mass, /r e — 2/(1 + A) 
is the mean electron molecular weight, and X is the H fraction. A 
small correction due to velocity reaching the sound speed rather 
than being null at the stellar surface has been neglected. 

Assuming the opacity for a given photon energy is a spatially 
independent constant times the density (as is the case, e.g., for 
Compton scattering or for bound-free absorption with a constant 
ionization coefficient), the optical depth is an integral over the pho- 
ton path, I (Pringle 1974), 

t(0)=to I d(/M (8) 


where r o = (<r c + cr lOD )noa, cr c is the Compton cross section (in- 
cluding its Klein-Nishina decline), cri 0n is a photoionization cross 
section per electron, n 0 is a fiducial electron density at r = a un- 
der the assumption of v = v m and a is the orbital separation. For a 
circular orbit and assuming the radiation originating at the compact 
object location, r is related to l by 


(r la) 1 = 1 + {l la) 2 - 2(1 /a) sin i cos </>, 


(9) 


where i is the binary inclination and </> is the orbital phase (defined 
as </> = 0, n at the superior and inferior conjunction, respectively), 
see fig. 1 in Zdziarski (2012). For r » R„ i.e., v ^ v<x,, 


n/2 + arcsin 5 

r(</>) - To -== — , S' = sin i cos </>, 

Vl - s 2 


( 10 ) 


which is is equivalent to equation (8) of Pringle (1974). 

In Compton scattering, a photon removed from one line of 
sight appears in a different one. Photons scattered away from di- 
rections with high t preferentially leave the system close to the 
directions with low r. Also, the photon distribution will be a func- 
tion of the angle with respect to the normal to the binary plane. It 
is possible to model this effect quantitatively in a more detailed 
treatment, e.g., using a Monte Carlo method. Such method was 
used by Hertz, Joss & Rappaport (1978), who, however, presented 
their results for i = 90° only (and for isotropic wind). However, 
the wind density distribution in Cyg X-3 is certainly not isotropic, 
but instead with a strong focusing towards the compact object 
(Friend & Castor 1982), and the wind structure is affected by ir- 
radiation from the X-ray source. Taking into account these compli- 
cations is beyond the scope of this work. We therefore make a sim- 
plifying assumption (as in Pringle 1974) that a scattering removes 
the photon from the line of sight and neglect photons scattered into 
it. This allows us to fit the observed phase-dependent fluxes by 


F(<f>)/Fo = exp [-r(0 - 0o )] • 


( 11 ) 


with r 0 , i, P, the offset phase, 0 O > and F 0 being the model free pa- 
rameters. We note that Fq = F max exp[r(7r)]. This requires the scat- 
tering optical depth to be <1 (though the absorption optical depth 
can be any). We define the fractional modulation depth as 

D = max — = 1 - exp(-Ar), Ar = r(0) - r(7r), (12) 

Fi nax 

where -F max and F m i n are the maximum and minimum modula- 
tion fluxes, respectively. For the model of equation (10), Ar ^ 
2roz'/ cos i. 

We then need to specify the parameters of Cyg X-3 appear- 
ing in this model. Since Cyg X-3 is a He star, X = 0 and /r e - 2. 
The stellar radius, R,, is closely related to the WR mass, see fig. 2 
of Schaerer & Maeder (1992). The separation, a follows from the 
total mass of the system. This, however, remains rather uncertain. 
Here, we choose, somewhat arbitrarily, a binary solution satisfying 
all the constraints of Vilhu et al. (2009), namely M, = 20Mq and 
the compact object mass of M c = 6Mq, for which a - 3.0 x 10 u 
cm. This M, corresponds to R, - 1 .Ox 10 1 1 cm (Schaerer & Maeder 
1992). This solution corresponds to an allowed range of inclina- 
tions of about i - (33°-63°), see section 6 of Vilhu et al. (2009). 
A typical value of the terminal wind velocity estimated for Cyg 
X-3 is Vo, - 1.7 x 10 8 cm s _1 (see a discussion and references in 
SZ08). We assume /? = 2, which was found to fit well the results of 
wind particle simulations for Cyg X-3 of Vilhu et al. (2009). Also, 
this value was used by Langer (1989) in his modelling of mass loss 
from WR stars. 

As an alternative, and purely phenomenological, model for the 
folded light curves, we use a sum of three Fourier harmonics (in the 
logarithmic space). We follow here exactly the treatment described 
in section 3.2 of Lachowicz et al. (2006). The modulation depth in 
this model is determined numerically. 
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Figure 3. (a) Folded/averaged light curves for the count rates for the energy 
range of 1.5-12 keV (blue, ASM), 14-50 keV (red, BAT) and 20^10 keV 
(green, INTEGRAL/ISGR1), using K = 40, 20 and 12 bins, respectively. For 
the BAT data, only observations with the exposure < 0.5 P/K have been se- 
lected. Here and elsewhere, two orbital ranges are shown for clarity of the 
display, (b) Comparison of the BAT folded light curve with the above expo- 
sure criterion and K = 20 with that without it (for K = 40). (c) Comparison 
of the ASM folded light curve calculated with our treatment described in 
Section 2.2 (blue triangles), i.e., using the barycentric-corrected data, the 
quadratic ephemeris of equation (3), the renormalization of the light curve 
to the running average, and the pre-averaging, with that obtained using a 
constant orbital period, unweighted averaging and no other corrections (red 
squares). 


2.3 Energy-dependent folded light curves 

We first consider folded/averaged light curves averaged over all the 
available ASM and BAT data. We have compared the light curves 
folded with the ephemerides of equations (2) and (3). We have 
found the differences are tiny, and we use the updated ephemeris 
of equation (3) hereafter. 

Fig. 3(a) shows the folded light curves for the total count rates 
of the ASM and BAT instruments (and also shows the 20^10 keV 
profile from INTEGRAL, see Section 3 below). We see the frac- 




Figure 4. Folded light curves using all of the data for (a) 1 .5-3 keV ASM, 
and (b) 14—20 keV BAT data, using K = 40. The red solid curves show fits 
to the data with the model of equation (8), (11). The green dashed curves 
show fits by the 3-harmonic model. 
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Figure 5. Comparison of the folded light curves in the hard (blue triangles), 
intermediate (green squares) and soft (red circles) spectral states for (a) the 
ASM 1.5-12 keV range, and (b) the BAT 14-50 keV range. 
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Figure 6. Comparison of the folded light curves in different energy channels 
for all the available ASM and BAT data. The blue, red and green light curves 
are for channels (a) 1 .5-3 keV, 3-5 keV, 5-12 keV, (b) 5-12 keV, 14-20 keV, 
50-75 keV, and (c) 24-35 keV, 75-100 keV, 100-150 keV, respectively. 


tional modulation depth, equation (12), is significantly larger in the 
1.5-12 keV range, D - 0.61, than the 14-50 keV one, D - 0.50. 
This is likely due to the bound-free absorption by the wind being 
stronger in the ASM band than in the BAT one. The BAT minimum 
is at 0/ 2n - 0.95, somewhat less than the ASM one at = 0.97. 
Both profiles have relatively complex and asymmetric shapes, with 
a flattening around the phase of 0.4 followed by a peak around 0.55- 
0.60. The ASM profile has two statistically significant dips, at the 
phases 0.40 and 0.575. For either of the dips, the difference with 
respect to the rate averaged over the two neighbouring bins is > 5 
standard deviations. Their presence indicates some complexity of 
the wind structure, probably due to focusing towards the compact 
object (Friend & Castor 1982). 

An issue for the BAT data concerns their relatively long typi- 
cal exposure of a single data point of ~ 10 3 s. This is comparable 
to the length of one orbital bin in the case of K = 20, which may 
cause a smoothing of sharp features in the folded light curves. In 


Fig. 3(a), we have thus used only measurements with the exposure 
< 0.5 P/K for the BAT profile. However, using only BAT data with 
short exposure times strongly reduces the number of used measure- 
ments. In the present case, accepting only exposures < 0.5 P/K for 
K = 20 reduces the number of available measurements by a factor 
of ~ 10 with respect to using all the available data. In Fig. 3(b), we 
compare the profiles with and without imposing this condition. We 
see that using the data without this selection causes certain smooth- 
ing of the profile, but the effect is rather minor. Thus, we do not use 
this selection for the BAT data hereafter. 

Fig. 3(c) illustrates how taking into account various effects can 
affect folded/averaged light curves, for the example of the ASM 
light curve. The red squares show the profile obtained with a con- 
stant P computed for the ASM data up to 2004 by Wen et al. (2006) 
(i.e., neglecting P), neglecting the barycentric correction, without 
prior renormalizing with respect to the running average and with- 
out pre-averaging within time bins. We see that that curve has the 
minimum at 4>l2n — 0.9 instead of 0.95 found for the treatment 
of Section 2.2. As we have checked, this shift is almost entirely 
due to neglecting the P of Cyg X-3. Taking into account barycen- 
tric time delays introduce only a very small effect. Furthermore, the 
red points have a much flatter and more scattered top part (similar 
to the profile shown in Wen et al. 2006). This difference is mostly 
due to renormalizing the light curve w/r to the running average. The 
profile without it is strongly dominated by the brightest states (as 
pointed out in Section 2.2). Then, fluctuations of those dominant 
states have a strong effect. On the other hand, our profile, shown by 
the blue triangles, is evenly averaged over all flux states. Finally, 
the pre-averaging introduces only a small effect for the ASM data, 
though it does noticeably increase the error bars. 

The 14-50 keV BAT profile shown in Fig. 3 can be com- 
pared with the folded light curve in the 12-50 keV band from the 
Gamma-ray Burst Monitor (GBM) on board Fermi, shown in fig. 
1 1 of Wilson-Hodge et al. (2012). Unlike our, rather smooth, BAT 
profile, the GBM curve has multiple sharp peaks. We note that the 
orbital profile for Cyg X-l shown in Wilson-Hodge et al. (2012) 
also has much stronger sharp features than the corresponding 15- 
50 keV BAT profile presented by Zdziarski et al. (2011). Possibly, 
the sharp GBM features are due to statistical fluctuations at some- 
what underestimated measurement uncertainties. 

In Fig. 4(a-b), we fit the lowest-energy channels of the ASM 
and BAT with the model of equations (8), (11). We see that neither 
profile is well fitted by this model. The profiles are significantly 
asymmetric, with a slower rise and a faster decline. (The scatter 
around the top of the 1.5-3 keV profile appears statistical.) Thus, 
the model does not provide statistically good fits, but we give the fit 
parameters for completeness: i = 33°±1°, 39°±1°, r 0 = 0.27±0.01, 
0.15 ± 0.01, 0 o /2zr = 0.02 ± 0.01, -0.03 ± 0.01, for the ASM and 
BAT, respectively (for the fixed /3 = 2, see Section 2.2). Given 
the statistically unsatisfactory fits, these parameters may not cor- 
respond to the actual wind of Cyg X-3. Also, the pairs of the pa- 
rameters (to, fi) and (;', fi) are each strongly anti-correlated. Still, 
our results indicate that, in spite of the strong X-ray orbital modu- 
lation, the inclination required to account for it does not need to be 
large, even for our spherically-symmetric wind model. The optical 
depths at the phases of 0, 0.25 and 0.5 are r = (1.5, 0.7, 0.5) and 
(1.0, 0.4, 0.3), respectively. Based on the fits, we can also calcu- 
late the corresponding mass loss rates using equation (7), the def- 
inition of r 0 , and the Cyg X-3 parameters adopted in Section 2.2. 
For the BAT and ASM best fits, r 0 corresponds to —M - (0.8- 
1.4) x 10~ 5 [trT/(o"c + cr ion )]Mg yr _1 (where cr T is the Thomson 
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Figure 7. The fractional modulation depth, D(E), for all the available ASM 
and BAT data (blue and black symbols) and the ISGRI data (magenta and 
cyan symbols) estimated using fits with the model of equations (8), (11) 
(blue and magenta) and with the 3-harmonic model (black and cyan). The 
dashed curve shows the dependence of D(E) for equation (7) in the case of 
the Klein-Nishina cross section, approximately normalized to 77(20 keV). 


cross section), which is of the order of the values usually estimated 
for Cyg X-3 (see SZ08 and references therein). 

F ig. 5(a-b) compares the folded light curves for the count rates 
of the ASM and BAT instruments for the hard, intermediate and 
soft states, using the criteria described in Section 2.1, see also Fig. 
2. We see that in the ASM range, the soft state modulation is sig- 
nificantly stronger than in the hard state, with the fractional modu- 
lation depth of Z) - 0.64, 0.57, respectively. A similar result for the 
ASM data (but using different state criteria) is shown in fig. 6 of 
Hjalmarsdotter et al. (2008). In the BAT range, the difference be- 
tween the states is smaller, D = 0.49 and 0.53 in the hard and soft 
state, respectively. The difference between the states may be due 
to the soft state taking place during periods of a higher wind mass 
loss rate. For the ASM and BAT, the ratio of the fitted values of To 
of our wind model between the soft and hard states is — 1 .20 and 
— 1 .08, respectively. The latter value may be close to the ratio of the 
mass loss rates (as the dominant opacity in the BAT range is likely 
to be Compton), whereas the former is likely to be also affected by 
changing the bound- free opacity due to the change of the ionizing 
continuum. 

We then compare the orbital modulation in different energy 
channels of the ASM and BAT, for all of the available data. We 
number the channels consecutively with the increasing energy, 
from 1 to 11. Fig. 6 shows some examples of comparison of the 
modulation profiles in different energy channels, and its panels are 
marked with the channel numbers being compared. We see that al- 
though the modulation depth generally decreases with the increas- 
ing energy, the strongest modulation occurs for the 3-5 keV range, 
as shown in Fig. 6(a). We have checked that the same behaviour 
takes place within the hard and soft states. This effect appears to be 
due to an additional, less or differently modulated, spectral com- 
ponent appearing in soft X-rays, which is probably re-emission in 
soft X-ray lines of the absorbed continuum by the stellar wind, 
see SZ08. The general decrease appears then to be due to the 
bound- free absorption cross section decreasing (except for ioniza- 
tion edges) and the Compton scattering cross section decreasing 
as well with the increasing energy. The profiles also become more 
symmetric with the increasing energy, indicating that the electron 
density is more symmetric (with respect to the plane perpendicular 
to the binary plane along the line joining the stars) than the ioniza- 




$/2n 
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Figure 8. Comparison of the folded light curves in different energy chan- 
nels selected by spectral state. The blue, red and green light curves are for 
channels (a) 24-35 keV, 50-75 keV, 75-100 keV in the hard state, (b) 24-35 
keV, 50-75 keV, 100-150 keV in the soft state, and (c) 24-35 keV, 35-50 
keV, 50-75 keV during periods of high-energy y - ray emission, respectively. 


tion structure. All the profiles have the minima within (f>!2n — 0.9- 
1, with no indication of a different modulation pattern up to 150 
keV. 

Fig. 7 shows the fractional modulation depth, D of equation 
(12), as a function of the photon energy, estimated using fits with 
the model of equations (8), (1 1) with a free and of the 3-harmonic 
model. Since the obtained statistical errors are very small, a mea- 
sure of the uncertainty is provided by the differences between the 
values for the two models. We show D only up to 75 keV, as at 
higher energies strong noise prevents trustworthy estimates. We see 
that indeed the orbital modulation is the strongest in the 3-5 keV 
range, and decreases towards both lower and higher energies. Fig. 
7 also shows D from the ISGRI data (Section 3 below), which is in 
a good agreement with the BAT results. 

The dashed curve in Fig. 7 shows the dependence predicted 
using equation (12) for Compton scattering alone, using the Klein- 
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Figure 9. Comparison of the folded light curves in the hard X-rays (red 
symbols) in the soft state in the 60-150 keV range with those at > 0.1 GeV 
(from FLC09, blue symbols). The constant background level (FLC09) has 
been subtracted from the LAT fluxes. The black curve shows the 1 00 keV 
modulation predicted (Z12) using the anisotropic Compton model fitted by 
Dubus et al. (2010) to the GeV data. All the profiles are normalized by the 
respective average values. 

Nishina cross section and Att = 0.75, for which the model de- 
pendence crosses the observed one around 20 keV. We see that the 
decline of D above 5 keV cannot be accounted for by the decline 
of the Klein-Nishina cross section with energy, and photoionization 
absorption has to contribute at <20 keV (in agreement with, e.g., 
SZ08). 

Figs. 8(a-b) show the modulation-profile dependence on en- 
ergy during the hard and soft state, respectively. We see that up to 
— 100-150 keV, at which energies the statistics become poor, the 
profiles have the minima around 0.9-1 phase, the same as at low 
energies and as for the entire data. An analogous figure for the in- 
termediate state looks similar and it is not shown. 

Fig. 8(c) shows the modulation for the data simultaneous with 
the observations of high-energy y-ray emission (see Section 2.1). 
We see the same modulation pattern as for other data, and, in par- 
ticular, neither a change of the positions of the maximum and min- 
imum nor a flattening of the profile with respect to the profiles at 
other states up to at least 75 keV. At higher energies, we also see 
no indication of the presence of another modulation pattern, but the 
data become noisy. 

The blue symbols in Fig. 9 show the orbital modulation ob- 
served at E > 0. 1 GeV by the Fermi LAT (FLC09). The black curve 
shows the modulation of the jet emission predicted at 100 keV us- 
ing the anisotropic Compton model (Dubus et al. 2010), from fig. 6 
in Z12. The 100 keV model is relatively similar to that at 0.2 GeV, 
though somewhat shifted to higher phases (Z12). These jet modu- 
lation profiles are shifted in phase with respect to those observed 
in X-rays, and have the maximum around the superior conjunction. 
We compare then the predicted modulation of the jet spectral com- 
ponent with that observed at 60-150 keV, red symbols. This profile 
has been obtained by averaging the BAT 100-150 keV profile with 
the INTEGRAL one for 60-100 keV (see Section 3 below), both in 
the soft state. Indeed, we see it is similar to the modulation pro- 


Table 1. The log of the INTEGRAL ISGRI observations. The total exposure 
is about 4.5 X 10 6 s. The first row includes all the observations before the 
start of the BAT monitoring. The low 14-50 keV state corresponds to the 
soft state as defined in Section 2. 


Start [MJD] 

End [MJD] 

Exposure [s] 

14-50 keV state 

52629.632 

53353.532 

2030100 

mixed 

53365.239 

53745.081 

32748 

high 

53750.259 

53881.889 

152052 

low 

53899.027 

53912.449 

4349 

high 

53912.449 

54059.913 

402936 

low 

54241.896 

54241.986 

6586 

low 

54438.871 

54444.899 

126269 

high 

54445.077 

54453.108 

168612 

low 

54574.516 

54609.314 

222279 

low 

54771.169 

54794.013 

341155 

low 

54824.233 

54973.407 

277455 

high 

55179.050 

55186.242 

69941 

high 

55341.028 

55341.047 

1796 

low 

55500.860 

55542.037 

651859 

high 


files observed in other states and very different from either the jet 
contribution predicted at 100 keV or observed at > 0.1 GeV. 

We note that Dubus et al. (2010) have determined the location 
of the y-ray emitting region at a distance ~ 2 a along the jet orig- 
inating at the compact object. Then, X-ray emission from the jet 
would be much less absorbed by the wind than that originating at 
the compact object. Furthermore, if the jet X-ray emission were 
stronger in the soft state then the modulation amplitude would then 
decrease. Instead, we see in Fig. 8 that the modulation depth ei- 
ther in the soft state or for the data simultaneous with the observed 
y-ray emission is basically the same as in the hard state, during 
which there is no detectable jet contribution. We note that the mod- 
ulation pattern of the jet 100 keV emission can be flatter than that 
shown in Fig. 9 due to the electron cooling slower than at > 0.1 
GeV, causing the 100-keV emission region along the jet to be more 
extended than that in high-energy y-rays. Still, a major jet contri- 
bution would flatten the observed total modulation profiles, which 
is not seen. Thus, based on Figs. 8-9 (also Section 3 below), we see 
no indication of a contribution from the low-energy tail of the y-ray 
emission. That emission has to have a low-energy cutoff above the 
hard X-ray range. Indeed, the jet models with the minimum Lorentz 
factors of the accelerated electrons of 1300-1500, which are shown 
by the blue curves in figs. 5(a-b) in Z12, do satisfy this constraint 
and contribute very little at 100 keV. On the other hand, their mod- 
els with the minimum Lorentz factors of 300 and 700 predict too 
much flux at 100 keV and appear to be ruled out. 


3 ORBITAL MODULATION IN INTEGRAL DATA 

We have studied INTEGRAL data from the ISGRI detector (a part of 
the IBIS telescope, Ubertini et al. 2003). We have used all the data 
public as of 2012 May 1, available at the INTEGRAL Science Data 
Centre (ISDC). We select only data with the off-axis angle < 5°, 
which corresponds to the fully coded part of the detector (the inner 
10° x 10° of the field of view). These data span MJD 52629-55542. 
INTEGRAL observations consist of pointings with a typical length 
of 2-3 ks. Our data set comprises of 1824 such pointings, for a total 
exposure of - 4.5 x 10 6 s. We consider only the 20-100 keV energy 
range since the data above 100 keV have relatively poor statistics. 
The data have been reduced using the Offline Scientific Analysis v. 
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Figure 10. Folded light curves for INTEGRAL/lSGRl observations, taking 
into account (a) all observations and (b) those corresponding to the soft 
state. The blue, red and green light curves are for energy ranges of 20-40 
keV, 40-60 keV, 60-100 keV, respectively. 

9.0 package provided by the ISDC (Courvoisier et al. 2003), with 
the pipeline parameters set to the default values. The created ISGRI 
light curves have a resolution of 100 s. Table 1 gives the log of the 
observations. 

We consider first all the observations, and then only those cor- 
responding to the soft-state intervals with the low BAT flux, see 
Section 2.1. We perform the analysis using the standard f tools 
software, in which the folded light curves are unweighted averages 
of the fluxes within individual bins. Our results are presented in Fig. 
10. We see no indication of a change of the profile shape towards 
high energies in either of the two cases considered, confirming our 
results in Section 2.3 above. We have also calculated the modula- 
tion profiles for only the data simultaneous with y-ray detections. 
We have obtained results entirely consistent with those for the soft 
state, but with somewhat worse statistics. 

In Fig. 3(a), we compare the 20^10 keV modulation profile 
for all the ISGRI data with the corresponding BAT profile. We see 
a very good agreement between them, which is consistent with the 
mostly overlapping energy coverage of the analyzed data. A 20- 
40 keV modulation profile for an early observation of Cyg X-3 by 
ISGRI is presented in Vilhu et al. (2003), whose results are similar 
to ours but with substantially higher statistical fluctuations. That 
paper also presents the modulation profile from the monitoring by 
the CGRO/BATSE during 1991-2000, which is also of almost the 
same shape as those of our ISGRI and BAT profiles. 

4 RXTE PCA/HEXTE PHASE-RESOLVED SPECTRA 

Here, we study the RXTE Proportional Counter Array (PCA) and 
High Energy X-ray Transient Experiment (HEXTE) data sets used 


Table 2. The log of the RXTE observations. The exposures are for the sums 
of all the used observations within each phase interval within a state. 


Spectral state 

Inferior Conjunction 

Superior 

Conjunction 


Exposure [s] 

Exposure [s] 


PCA 

HEXTE 

PCA 

HEXTE 

i 

9120 

9120 

23408 

19549 

2 

24048 

14611 

42976 

25952 

3 

20912 

23049 

23008 

9337 

4 

3744 

4645 

2960 

3805 

5 

1120 

5465 

12688 

11672 


by SZM08, who obtained average X-ray spectra in five spectral 
states. The log of the observations is given in Hjalmarsdotter et al. 
(2009). We split each average data set into two subsets correspond- 
ing to the superior and inferior conjunctions, which is done by di- 
viding the orbital phases into two parts, 0/2;r = 0.28-0.78 and 
0.78-0.28. These boundaries have been determined using the or- 
bital template of van der Klis & Bonnet-Bidaud (1989) with the 
condition that the orbital period is divided into two parts of equal 
length, with the template flux greater or lower than that correspond- 
ing to the dividing phases. 

We have performed standard RXTE PCA and HEXTE data re- 
duction using ftools. The resulting exposures for the two subsets 
of the data of SZM08 are given in Table 2. We see that whereas the 
states 1-3 have exposures > 10 4 s for each subset for either PCA 
or HEXTE, the states 4-5 have much worse coverage, especially 
for the PCA, with the exposures as short as - 1000 s and - 3000 
s. Therefore, we present only the spectral shapes for the states 4- 
5, whereas we study in detail the differences between the subsets 
around the superior and inferior conjunctions for the states 1-3. 

In our approach, we study the difference between the spectra 
around the inferior and superior conjunctions assuming it is due to 
an additional optical depth, with scattering and absorption charac- 
terized by a constant ionization coefficient, see equation (B3) of 
Zdziarski (2012). We first fit the spectra around the inferior con- 
junction. Finding a physical model of unabsorbed spectra of Cyg 
X-3 has proven, in this and previous studies, to be difficult. For ex- 
ample, Hjalmarsdotter et al. (2008, 2009) found that fitting data at 
> 3 keV do not allow an unambiguous determination of the intrinsic 
spectra, with the range of possible models strongly varying at low 
energies. On the other hand, Cyg X-3 data from BeppoSAX extend 
to < 1 keV, and are better suited for such a determination. Those 
data have been studied by SZ08, who used a relatively complex 
model of wind absorption and emission, modifying spectra from 
hybrid plasma and Compton reflection (Poutanen & Coppi 1998; 
Coppi 1999; Gierlinski et al. 1999). 

Here, we fit the data from PCA and HEXTE, which are for 
E > 3 keV only, and do not allow unambiguous determination of 
the low-energy intrinsic spectra. However, our objective is to study 
physical mechanisms of the orbital modulation, and determining 
the actual intrinsic spectral distribution is beyond the scope of this 
work. We thus fit the data with a model that has been shown to pro- 
vide relatively good fits to a large set of data from RXTE and INTE- 
GRAL (Vilhu et al. 2003; Hjalmarsdotter et al. 2008, 2009, SZM08, 
Corbel et al. 2012), but which still should be considered as phe- 
nomenological. In particular, any inferences about the unabsorbed 
spectra based on that model are very uncertain. The model is based 
on emission of hybrid plasma and reflection (including the Fe K 
line), similar to the model used by SZ08. However, the complex 
absorption by the wind is treated only phenomenologically, as a 
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Figure 11. The panels (a-c) show the spectra for the states 1 -3, respectively, 
in the classification of SZM08. The upper (blue) and lower (red) spectra cor- 
respond to the phase ranges around the inferior and superior conjunction, 
respectively. The crosses show the unfolded PCA+HEXTE spectra. The in- 
ferior spectra are modelled by emission of hybrid plasma and reflection/Fe 
K line absorbed by a dual absorber. The models of the superior spectra are 
the best-fit models of the inferior spectra absorbed/scattered in an additional 
medium, corresponding to the difference in the average column density be- 
tween the two phase ranges. 


product of full and partial neutral absorbers. The Fe abundance of 
the absorber and reflector is assumed to be free. However, given the 
neglect of photoionization and the complex spatial structure of the 
wind, the fitted value should also be considered as phenomenolog- 
ical only. In spectral fitting, we allow a free normalization of each 
HEXTE spectrum with respect to the corresponding PCA one, and 



Figure 12. The spectra of the states 4-5, shown in the panels (a-b), respec- 
tively. The upper (blue) and lower (red) spectra correspond to the phases 
around the inferior conjunction and superior conjunction, respectively. The 
crosses represent the unfolded PCA+E1EXTE spectra. The models (not 
shown) are given by emission of hybrid plasma and reflection/Fe K line 
absorbed by a dual absorber. 


we add a 1 per cent systematic error to each of the count spectra. 
The PCA and HEXTE data are fitted in the energy ranges of about 
3-25 keV and 20-40 keV, respectively (as we have found the sta- 
tistical quality of the HEXTE spectra above 40 keV is relatively 
poor). The resulting fits to the spectra in the states 1-3 are sta- 
tistically satisfactory, X 2 1 V< Z * ■ The fitted parameters are consistent 
with those presented in Vilhu et al. (2003) and Hjalmarsdotter et al. 
(2008, 2009). 

Then, we take the best-fit model to an inferior-conjunction 
spectrum, and fit the corresponding superior-conjunction one with 
the inferior-conjunction spectral model but attenuated by absorp- 
tion and scattering in an additional medium due to lines of sight 
through the wind around the superior conjunction being longer 
than those around the inferior conjunction. We have found that the 
fractional modulation depth averaged over our phase intervals is 
D = 0.4 but it strongly decreases at E> 15 keV. We have first 
tested whether this behaviour can be due to the Klein-Nishina de- 
cline of the Compton cross section. However, that decline is much 
too slow to account for the observed energy-dependent modulation 
depth, which confirms our results in Section 2.3 above. Clearly, 
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Table 3. The fitted parameters of the additional absorbing medium related 
to the superior conjunction for the states 1-3 (which have the exposures 
sufficient for averaging). See Section 4 for details. 


State 

Afuii 
1 0 23 cm 2 

^part 

10 23 cm 2 

fcaw 

logio £ 
erg cm s -1 

x 2 /y 

t 

2.0 ±0.2 

6.8 ±0.5 

0.36 ±0.01 

1.6 ±0.4 

40/53 

2 

2.1 ± 1.4 

10.8 ±8.7 

0.31 ±0.08 

3.3 ±0.4 

30/55 

3 

3.5 ±0.3 

3.0 ±0.6 

0.44 ± 0.07 

1.8 ±0.4 

55/56 


photoionization, known to be present in Cyg X-3 (e.g., SZ08), also 
contributes. Its cross section generally declines fast with increas- 
ing energy (apart from ionization edges). Accordingly, we have 
found that good fits to the data can be obtained including absorp- 
tion by an ionized medium. Specifically, our model for the addi- 
tional superior-conjunction attenuation consists of partial covering 
by a partially ionized medium and full covering by a fully ionized 
medium. Such a model may approximate shadowing of the stellar 
wind by the X-ray source, and it has been invoked to explain the 
IR modulation and spectral lines (van Kerkwijk 1993). The attenu- 
ation is then by the factor, 

A(E) = exp(-(V ful |<T c ) { 1 - / cov + / cov exp [-A part (cr c + 

^"ion )]}.(13) 

where Nf^a and N par t are the electron column densities of the fully 
and partially covering medium, and f cow is the covering fraction. 
We use the ionized absorber of Reeves et al. (2008), based on 
xstar (Bautista & Kallman 2001). It is characterized by the col- 
umn density and ionization parameter, £ = L x /(nr 2 ), where n is the 
medium density and r is the distance from the ionizing source to the 
medium; see Reeves et al. (2008) for further details. Our model is 
clearly very simplified. We have neglected spatial variations of the 
ionization coefficient. Also, we have assumed that scattered pho- 
tons are entirely removed from the line of sight, whereas some pho- 
tons are also scattered into the line of sight from other directions. 
The fit results are summarized in Table 3, and the unfolded spectra 
and the models are shown in Fig. 11. 

We see that the total electron column density is N ~ 10 24 cm 2 , 
which corresponds to the difference in the average Thomson optical 
depth between the two phase intervals of Att ~ 1 . This is in agree- 
ment with Ar ~ 1 found from fitting the orbital modulation profiles 
(Section 2.3). The fitted values of £ ~ 10 2 ~ 3 erg cm s _1 correspond 
to moderate ionization. We can notice that the models do not fit the 
smeared edges clearly seen around ~ 10 keV, which are apparently 
due to absorption by H-like Fe in a very strongly (but not fully) 
ionized phase of the wind, which we do not model here. Using the 
definition of £, the characteristic distance of the ionized medium is 
r ~ L x /(Nt;), which, for the fitted values of £ and N is r ~ 10 11-12 
cm. This is of the order of the orbital separation, a, as expected. 
Thus, the fitted model is consistent with bound- free absorption tak- 
ing place in the stellar wind. We note that the fitted parameters do 
not form monotonic sequences with the increasing state number. 
This appears to be due to the limited observation exposures. Then, 
the two sub-states of a given state are not true averages, but are 
affected by variability taking place during the used observations. 

In Fig. 12, we show the inferior and superior spectra for the 
states 4-5. We can see that due to the short exposures, the intrin- 
sic variability is apparently comparable to the orbital modulation, 
especially for the spectra of the state 4, which have the superior- 
phase fluxes at > 20 keV higher than those of the inferior phases, 
contrary to the expectations for long-term averages. 


5 CONCLUSIONS 

We have studied orbital modulation of X-rays from Cyg X-3. Us- 
ing the data from the detectors ASM, PCA and HEXTE of RXTE, 
BAT of Swift, and ISGRI of INTEGRAL, we have found that the 
modulation depth at E > 5 keV decreases with the photon energy. 
The decrease is too fast to be accounted for by the Klein-Nishina 
decrease of the Compton cross section with energy. As we confirm 
by fitting the PCA/HEXTE data, the energy dependence requires 
the presence of a moderately ionized absorber, with the ionization 
coefficient of L x /{nr 2 ) ~ 10 2 ~ 3 erg cm s _1 . We have also found 
a decrease of the depth of the modulation at E < 3 keV, which is 
probably due to the re-emission of the absorbed continuum in soft 
X-ray lines by the wind. 

Given that we use the data accumulated over long time, we 
determine the energy-dependent folded and averaged light curves 
with high accuracy. The modulation profiles are not consistent with 
a spherically symmetric wind. The measured detailed shapes of the 
orbital modulation and their photon-energy dependence can serve 
for determination of the wind structure in future work. Also, the 
energy dependence of the modulation depth at low energies can be 
further studied using spectra from detectors sensitive at soft X-rays. 

Assuming that the modulation is from wind absorption, the 
phase of the flux minimum corresponds to the maximum optical 
depth through the wind, which is then close to the superior conjunc- 
tion. On the other hand, high-energy y-rays detected from Cyg X-3 
by Fermi are found to have an orbital modulation shifted by about 
a half of the period with respect to that of the X-rays (FLC09). 
This pattern has been interpreted as due to anisotropy of scattering 
of stellar photons by relativistic electrons in the jet by Dubus et al. 
(2010). As calculated by Z12, the jet emission at hard X-rays has 
the modulation pattern almost the same as that in the GeV range. 
Thus, the jet X-ray contribution can potentially be measured by de- 
termining the orbital-modulation profiles at hard X-rays. 

We have performed such study, and have found no presence of 
an X-ray modulation pattern different from the standard one, i.e., 
with the minimum around the superior conjunction, up to ~100 
keV. This implies that an X-ray jet contribution is weak up to at 
least that energy and the jet component has to have a low-energy 
cutoff. This is consistent with the X-ray spectral measurements 
taken during periods of y-ray emission, which do not show a hint 
of a jet component up to comparable energies (Z12; Corbel et al. 
2012). This constraint is satisfied by two among the models pre- 
sented in Z12, those with the minimum Lorentz factors of the ac- 
celerated electrons of 1300 and 1500. 
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